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The origin of the breakdown scale in an effective field theory treatment of nuclear 
forces is investigated. Different organizational schemes used in the nonperturbative 
calculation of nucleon-nucleon scattering seem to break down for different physical 
reasons. It is argued that the same physics is actually responsible, namely strong 
two-pion exchange as also observed in the pion scalar form factor. The prospect 
of extending the applicability of the effective field theory beyond this breakdown 
point is discussed. 



1 Introduction 

Strong interactions at low energies can be described by an effective field tin 
ory (EFT) of pions and nucleons with the nucleons treated nonrelativistically. 
Chiral symmetry prescribes how pions interact. All other interactions are pa- 
rameterized by the most general local Lagrangian consistent with the symme- 
tries of quantum chromodynamics (QCD). 

The coefficients of these general terms can either be predicted directly 
from QCD, which may be possible in future lattice calculations, or fit from 
existing data. I will focus on nucleon-nucleon (NN) scattering in the following. 
The data here is quite good and should lead to an excellent determination 
of these constants. However, since nuclear forces produce bound states, a 
nonperturbativa treatment of the EFT is required, which adds a new twist to 
this procedure. El 

Any EFT will break down at a scale associated with the underlying physics 
not taken into account in the original Lagrangian. For nuclear EFTs, an op- 
timist may believe this will occur when the center of mass momentum of the 
nucleons excite either the exchange of the next heaviest meson after the pion, 
p ~ rup, or the nucleon-A mass difference p = a/ M{Ma — M) ~ 525 MeV. 
In practice, the breakdown scale seems to be much lower, p 300 MeV. This 
could be due to anything from quark substructure to subtleties in the fitting. 

Furthermore, two different organizational schemes for calculations in the 
EFT seem to imply different physical reasons for the breakdown. This would 
be unfortunate, since alternative ways to organize the corrections should have 
no effect on the physics. By using some models, I will show what actually does 
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Figure 1: The analytic structure of an observable for the two- Yukawa model of Eq. (| 
the complex momentum plane. 



set this scale and then use this to make connections to the real world. 

Extending the applicability of the EFT to the scale associated with the A 
seems like a small achievement, but could make calculations of nuclear matter 
converge much faster. In that case, the momentum of interest is on the order 
of the Fermi momentum at saturation pf ^ 280 MeV, which currently is too 
close to the breakdown scale for an EFT analysis to be useful. 

2 Removing vs. Weakening the Pion 

Assume, for now, that interactions between two nucleons are exactly given 
by potential exchange of both a long- and short-ranged meson with masses 
TOtt — 140 MeV and rup = 770 MeV respectively. 
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The pion coupling is taken to be a-^r = ffA™^/(167r/^) ~ 0.075, as in the real 
world, and the rho coupling ap = 1.05 is tuned to give a large scattering length 
a = —23.4 fm, as observed in data for the partial wave of NN scattering. 
The following discussion does not depend on the value of the potential at the 
origin, so I will take Ksxact to be zero there for simplicity. 

This potential generates observables with analytic structure as shown in 
Fig. 1^, specified by a cut at both m-^jl and mp/2. The existence of these cuts 
indicates that a low-momentum expansion of an observable will not converge 
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for momenta above rfiT^jl. An EFT can only go beyond this point if it either 
exactly reproduces or sufficiently approximates the analytic structure of the 
pion. 

This can be illustrated by focusing on the most general Lagrangian without 
pions (I will restrict myself to the ^5*0 channel below) 

N ~ CoiN^Nf + [{N^VNY + h.c] +... , (2) 

which corresponds to a nonrelativistic potential between two nucleons of 

Veft{p,p') ^Co + {p" + p'^) + . . . , (3) 

where p and p' are the center of mass momentum of the incoming and outgoing 
nucleons. Since these interactions depend on powers of momentum, a naive 
calculation of the NN scattering phase shifts from CcS leads to ultraviolet 
divergences. This just reflects our ignorance of the high-energy behavior of 
this theory. 

Two possible procedures to regulate these divergences are: i) adding a 
cutoff to all integrals to eliminate or suppress the high momentum states and ii) 
using dimensional regularization, which amounts to throwing out any integrals 
that diverge like powers of momentum. The second procedure may seem a bit 
counterintuitive at first, but is simpler to implement since it does not break any 
symmetries of the underlying theory and therefore does not require additional 
counterterms. Perturbative calculations in chiral perturbation theory have 
confirmed the equivalence of these regularization procedures. 

Nonperturbative calculations, on the other hand, show that a naive im- 
plementation of dimensional regularization is incomplete. A generalization to 
account for integrals that linearly diverge is required. This regularization is 
called power divergence subtraction^ (PDS). 

After regulating the theory, the EFT can be organized to give predictions 
order-by-order in powers of an expansion parameter. Since the low-energy con- 
stants C2n are not determined by the symmetries of QCD, they must be fit to 
data. In the process, they will typically pick up the next scale A of physics 
not explicitly accounted for in the EFT, leading to C2„p^" ~ (47r/A^)(p/A)^". 
Therefore interactions with more derivatives are suppressed by additional pow- 
ers of p/A. Power counting in this parameter allows Eq. to be truncated 
to the first few terms, with more accurate results obtained by keeping more 
terms in the momentum expansion. 

There are two prevalent ways of implementing power counting in nonper- 
turbative effective field theories. The first is to count powers of A in the poten- 
tial Veft and then solve the Lippmann-Schwinger equation. Such a procedure 
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sums an infinite number of Feynman diagrams to generate the nonperturbative 
amplitude A, related to the phase shift S by 

. ■ (4) 

p cot — ip 

This method of power counting is called A-countingi and is in the spirit of t 
Hamiltonian theory of Wilson, first applied to nuclear EFTs by Weinberg 
An alternative is to expand A in powers of Q = (p, to^, 1/a), motivated by 
chiral perturbation theory, which expands in p and m^. The nonperturbative 
nature of the interactions identifies 1/a as a small parameter as well, justifying 
its inclusiDH in the expansion. This method of power counting is called Q- 
counting. Bb 

For a theory without pious, such as Eq. (^), the underlying scale A is just 
TOtt. Both power counting procedures work as long as momenta are less than 
this scale, reproducing the effective range expansion 

pcot,5 = -- + -p2^r„^. (5) 

n=0 

The effective ranges, which are determined from NN scattering data, i also 
exhibit this scale r„ ^ 1/A ^ l/mjr. To extend the EFT to higher momenta, 
pions need to be explicitly added to the effective Lagrangian. This hopefully 
will bring A up to the QCD scale of rUp or possibly even 1 GeV. 

The inclusion of pions nonperturbatively is different in the two power 
counting procedures. Every iteration of one-pion exchange contributes to the 
pion cut in Fig. |l|. In general, n iterations are proportional to^ 
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where Ann — 300 MeV is just another way of expressing the strength of the 
interaction via the relation = to^/(MAnn)- A-counting includes this long- 
range part of one-pion exchange fully at leading order, removing the effect 
of the pion cut. Q-counting pushes one-pion exchange to subleading order 
and iterated one-pion exchange to even higher orders. Therefore, each order 
in the Q-counting expansion weakens the pion cut by an additional power of 
m,r/ANN, as seen in Eq. (^). The topic of the next section is whether this 
simpler implementation of power counting, also called perturbative pions, can 
be as successful as A-counting. 

"^This can be seen by studying the exact solution of an exponential potential, i which has the 
same cut structure as a Yukawa or by looking at the analytic eejpressions of higher Yukawa 
loops, which exhibit the cut structure through Polylogarithms. u 
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Figure 2: The error in pcot<5 for a pion coupling three times stronger than the real world. 
The left plot is from a normal fit to the low-energy constants both without (solid) and with 
(dashed) explicit pions. The right plot uses the modified effective range expansion to fit. 



3 New Fitting Procedures 

The NN EFT is expected to work best at low-momentum. So a natural pro- 
cedure to fix the low-energy constants C2n is to calculate an observable, like 
pcot(5, and fit near threshold. However, there are some complications originat- 
ing from the nonperturbative nature of the problem. These must be cleared 
up before being able to make a statement about the breakdown scale. 

3.1 A-counting and the Modified Effective Range Expansion 

Most calculations that implenient A-counting require a numerical treatment 
(although there are exceptionsB). A nice way to show the jSy^jtematic improve- 
ment expected of an EFT in this case is to use error plotsH'Qfl as in Fig. ^. As 
more and more constants are added to the effective potential Eq. (^, the error 
in p coil 5 improves systematically by powers of p^, reflected in the improving 
slopes. The point at which all errors are of the same order is where the theory 
is expected to break down. El The solid lines of the left plot of Fig. |2| show the 
case without explicit pions. The breakdown point, or radius of convergence, of 
the EFT is around as expected from the discussion in Section ^ 

Adding pions to the short-ranged Lagrangian Eq. as dictated by chiral 
symmetry should improve the radius of convergence. E3 However, one needs 
to be careful when implementing this with cutoff regularization. If the pion 
coupling were three times stronger, for example, adding it to the Lagrangian 
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Figure 3: Symbolic definition of quantities in the modified effective range function in Eq. 
which sums all the pion contributions. 



would not improve the breakdown scale. B This is illustrated by the dashed 
lines in the left plot of Fig. ^. Chiral symmetry helps in that the coupling 
is small, but an EFT should work regardless of the strength of the potential 
being added. What went wrong? 

Note that the effective ranges r„ of Eq. (||), which are used to fix the 
values of the C^n , exhibit the scale to^ . If this scale were completely removed 
from the fitting procedure, then the constants should take on the proper scale 
rup. One way to do this is to construct a function similar to pcotS which 
is analytic all the way to the rho cut in Fig. ^. In proton-proton scattering 
analyses, such a procedure is used to remove the threshold cuts coming from 
photon exchange. Called the modified effective range expansion, til it can be 
used here to remove the effect of pion exchange by defining 
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All possible insertions of one-pion exchange are taken into account in this 
formula, as seen pictorially in Fig. ||. For the model potential Eq. (|l|), this 
function is analytic up to mp/2 and has modified effective ranges rn ~ 1/A ^ 
l/nip. Fitting the C2n's using the function K{p) leads to much better results, 
shown in the right plot of Fig. |[ The radius of convergence is now on the 
order of rUp as expected. 

Using a cutoff to regulate the theory mixes loiig- and short-distance physics, 
which is what led to the premature breakdown. El The modified effective range 
expansion removes the long-distance physics from the observable, allowing a 
clean fit of the low-energy constants to the short-distance physics. In princi- 
ple, this can also be achieved by the addition of more counterterms, but the 
modified effective range expansion is a simpler alternative. Luckily, one-pion 
exchange in the real world is weak enough that an improper treatment of this 
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Figure 4: Comparison of the two power counting procedures for the two- Yukawa model 
Eq. using a fit at low momentum. 

malady does not make a difference, i However, this might not be the case for 
two-pion exchange as wiU be discussed in Section H. 



3.2 Q-counting and a Global Fit 

A low-momentum fit to data using Q-counting and perturbative pions also 
breaks down much earlier than expected. □ In Fig. ^, the leading order (LO) 
and next-to- leading order (NLO) calculations using Q-counting are compared 
to those using A-counting. The A-counting result has a radius of convergence 
around nip after the low-energy constants are matched to the modified effective 
range expansion. 

Implementing the modified effective range expansion is inconsistent with 
the idea of perturbative pions, and so the Q-counting result breaks down 
around m,r. However, it is important that the low-energy constants only pick 
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up the (unknown) short-distance physics, and matching to the modified effec- 
tive range expansion is the mpst convenient way to do this. Using PDS, this 
would lead to the expressions E£l 

which only depend on the short-distance physics.^ One important modification 
to this would be to match the pole of the amplitude A exactly l3 

+ 1/a 

Otherwise, terms attempting to correct the position of the pole would enter as 
double poles at higher orders messing up the power counting. This modifica- 
tion requires perturbative corrections to Co that compensate the pion contri- 
butions at each order. However, this still leaves the other C2n's dependent on 
the modified effective ranges r„, which cannot be obtained using perturbative 
pions. 

A method must therefore be devised to carry out the matching of Eq. (H). 
It is reasonable to assume that when the modified effective ranges are set to 
the correct values, the error of the EFT will be less than when the ranges are 
at some random values. This can be checked by using a model. Replacing the 
Yukawas in Eq. (|^) by delta-shells 

- > -gi Sir , (10) 



an analytic expression for the exact answer can be obtained. Expanding the 
exact amplitude in powers of Q, it can be compared to the blind EFT fit of 
the low-energy constants to data using Q-counting. The results coincide best 
when the modified effective ranges are varied to globally minimize the for 
data in a representative momentum range like [to^ , mp] . 

The results are shown up to next-to-next-to-leading order (NNLO) in 
Fig. ||. They suggest that for Q-counting, a constrained global fit of the EFT 
to data to fix the low-energy constants is better than a low-momentum fit. The 
constraint is that the modified effective ranges, not the low-energy constants 
themselves, are varied to achieve the best fit. 

One thing to notice about Fig. || is that the errors do not show the system- 
atic improvement in slope characteristic of the earlier error plots. Although the 

''Although the renormalization scale fi appears in these low-energy constants, all observables 
are independent of it. 
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Figure 5: Q-counting expansion of the two delta-shell model. The exact expression (solid) 
coincides nicely with the EFT (dashed) when a global fit is minimized with respect to the 
modified effective ranges. The dashed line at LO is on top of the solid line. 



scattering length is fit exactly to reproduce the pole in the amplitude Eq. 
the other terms in the effective range expansion have corrections at each order 
in the Q-counting. A comparison of Fig. ^ with the right plot of Fig. || shows 
that accuracy is sacrificed near threshold to obtain the global fit needed for this 
power counting. In addition, the ability to read the EFT radius of convergence 
from the error plot is gone. 

Another way to estimate the breakdown scale is therefore needed. The 
calculated amplitude A starts at order and has a contribution An at each 
order Q" in the expansion. Up to this piQint, I have extracted the Q-counting 
phase shift by expanding the expression^ 
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but I could also have kept the full expression Eq. (pT|). These two expressions 
are equivalent up to terms that are higher order in the Q expansion. Those 
extra terms, though, are an estimate of the corrections to the actual result 
and depend on the radius of convergence. The breakdown for Q-counting can 
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Figure 6; Determination of the Q-counting breakdown scale for the two- Yukawa model. 



be identified as the point at which the two expressions Eqs. (^Tj) and (12) 
diverge from each other. Doing this exercise for the two- Yukawa model results 
in Fig. ^. There it is easy to see the actual breakdown is somewhere around 
vip as expected. 



4 Data and the Role of Two-Pion Exchange 

Now that the two power counting methods are well understood and give the 
proper radius of convergence for a model, it is appropriate to apply them to 
actual NN-scattering data. Using the Nijmegen partial wave analysis, Q the 
radius of convergence for each of the two different power counting schemes 
is shown in Fig. ^ The left plot shows that, for A-counting with one-pion 
exchange included as in the previous section, the breakdown is around p ^ 
300 MeV. The right plot shows a similar scale for Q-counting with perturbative 
pions. 

This scale is less than expected and seems to be set by physics other than 
the p or A. Early calculations in Q-counting predicted a breakdown at Ann, 
which would coincide with the observed scale. However, recall that in the 
previous section I showed that Q-counting broke down around vrip in the two- 
Yukawa model. That calculation had identical content with respect to one-pion 
exchange, indicating that the 300 MeV comes from different physics. A closer 
study of A-counting determines that two-pion exchange should be included at 
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Figure 7: Results for ^So NN-scattering data. The breakdown scale is determined by the 
left plot for A-counting and the right plot for Q-counting. 

NLO.0 

The fact that the actual data behaves so much differently from the two- 
Yukawa model indicates the model was not sophisticated enough. Furthermore, 
no reasonable choices for the parameters of this model can give as large an 
effective range rg = 2.63 fm as seen in the ^5*0 channel. □ Until now I have 
assumed that and rup are the underlying scales, but nuclear physics wisdom 
tells us there are actually two short-range Yukawas 

g-m„r g-m^r ^-nipV 

V{r) ^ -a^r hap , (13) 

r r r 

given by an attractive scalar (of mass between 500 and 600 MeV) and repulsive 
vector particle.^ Taking the short-distance couplings to be much stronger than 
the pion (acr,ap) = (7[-|14.65), the modified effective ranges r„ of the data 
can also be-jreproduced.El This three- Yukawa model is reminiscent of the Bonn 
potential, llil which is known to model the data well. 

Applying the NN EFT to observables produced from Eq. ( [1^ ) indeed leads 
to a breakdown around ma/2 ~ 300 MeV. I can verify the breakdown is not set 
by the two-pion threshold 2m^ within this model by varying the pion and scalar 
masses. The (t(600) particle is a common parameterization in nuclear physics 
for two-pion exchange contributions. The large scalar coupling indicates that 
including such effects through chiral symmetry might require the modified 
effective range expansion. Such an analysis could improve the EFT radius of 
convergence.^ 

'^The vector particle is normally taken to be the lu meson, but for continuity with the earlier 
discussion, I will call it the p meson. 

''A preltoiinary calculation, inspired by discussions at this workshop, seems to indicate this 
is true.li^l 
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It might seem odd that two-pion effects could be strong, since they should 
be weak according to chiral symmetry. 13 However, there are other instances 
where two-pion effects have proven to be strong. The pion scalar form factor, 
defined by 

{7T''{p)Tr\p')\m{uu + dd) |0) = <5"V,r(s), s = (p + p')\ (14) 

has been shown to have poor convergence in chiral perturbation theory, 
breaking down just above threshold at ^/s ~ 400 MeV. This comes about 
from infrared effects that require the summation of (Ins)" terms to all orders. 
Carrying out this summation does lead to improved results and the generation 
of a "resonance" -like shape in the dispersion analysis just like a scalar resonance 
around 600 MeV. tj The same phenomenon could occur in NN scattering and 
requires further study. 

5 Red Herrings 

The models developed to understand the breakdown scale better can also be 
used to shed light on other convergence issues. One feature of the three- 
Yukawa model is that it can reproduce the somewhat large effective range in 
NN scattering, which is ro ~ Any reasonable choices of parameters for 
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the two- Yukawa model can only give ro around Therefore, it is possible 

that the lower radius of convergence of the three- Yukawa model might have 
nothing to do with two-pion effects, but come from some type of correlation 
with the effective range rg instead. Choosing different parameters in Eq. (p^, 
I can keep a, rg, and the mass difference rUp — fixed while changing 
from 500 MeV to 1500 MeV. The results are shown in Fig. ||, confirming that 
the radius of convergence is set by vfiajl alone. 

Recently, the convergence of Q-counting has been challenged within the 
context of a low-energy theorem. El The point can be made by defining the 
observable 

5(/)=pcot^- (^-i + irop2^ , (15) 
and taking -p — |l/a|. Then Q-counting has two types of corrections: 
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It was shown that the Q-counting prediction for S{p'^) at NLO is very far from 
the data. 

This can also be seen in a model. For an analytic understanding, I will 
use the two delta-shell potential of Section ^. Expanding in Q, 
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one finds the NNLO term is larger than the NLO term. Only after a few more 
orders does the series start to converge. 

Taking the exact answer, I can also expand in Q to see what is happening 
analytically. These expressions are messy and can be simplified by observing 
that even though m^ra = 0(1) in Q-counting, it is actually of order 10 nu- 
merically. Expanding 3(1/0^) in (m^a)~^ then should simplify the expression 
immensely without blurring its meaning, leading to 
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(19) 



Terms in the second bracket and higher are suppressed by (rriT^a) ^ and can be 
neglected. Focusing on the first bracket of terms, it is clear that the expansion 
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is in TOtt/Ann as expected for Q-counting. However, the leading term (which 
happens to be NLO for this particular observable) has a numerical factor of 
i, which suppresses it down to the order of the NNLO term. Therefore, the Q 
expansion of S{p^) works, but is slowly converging due to some inconvenient 
numerical factors as well as the size of the expansion parameter, Eq. ([l6|). 



6 Summary 

Nonperturbative effective field theories are subtle to implement. Techniques 
used successfully in chiral perturbation theory, such as naive dimensional reg- 
ularization and low-momentum fitting, do not produce the expected results in 
nonperturbative calculations. After a careful study using solvablejpaadels, new 
methods have been developed to deal with these complications. Hu'tHl 

Pions enter the calculation differently depending on whether power count- 
ing is done in the QCD scale (A-counting) or the chiral scale (Q-counting) . 
In both cases a proper determination of the low-energy constants from the 
short-distance physics requires more than just a low-momentum fit. Utilizing 
the modified effective range expansion to remove the effect of the long-distance 
pions from the procesa.of fitting, the A-counting radius of convergence moves 
beyond the pion scale. 13 A constrained global fit minimized with respect to t 
modified effective ranges is required to achieve the same effect in Q-counting. 

Applying the new fitting procedures to NN scattering data results in a 
breakdown of about 300 MeV. Models can be used to confirm this breakdown 
is associated with the same physics for both power counting procedures. Using 
a three- Yukawa model, it was argued that including two-pion exchange will 
extend the EFT radius of convergence. It is not unreasonable to assume t 
effect is strong since a similar situation occurs in the pion scalar form factor. li 

Finally, the convergence of Q-counting was investigated using an exactly 
solvable model. It was shown that some observables could have inconvenient 
numerical factors, such as occur in the low-energy theorem, □ which suppress 
the leading orders and mask the behavior of the expansion. However, analytical 
results show that this expansion does converge, albeit slowly. 
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